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We consider several processes operating during the late stages of planet 
formation that can affect observed orbital elements. Disk-planet interactions, 
tidal interactions with the central star, long term orbital instability and the 
Kozai mechanism are discussed. 



I. The formation of protoplanets 

Observations. It is generally accepted that the planets in our so- 
lar system were formed in a flattened gaseous nebula centered around 
the Sun. In typical star-forming molecular clouds, dense cores are ob- 
served to have specific angular momentum larger than 6 x 10 20 cm 2 s _1 
(Goodman et al. 1993) such that their collapse leads to rotationally 
supported disks analogous to the primordial solar nebula (Terebey et 
al. 1984). Between 25 to 75% of young stellar objects (YSOs) in the 
Orion nebula appear to have disks (Prosser et al. 1994; McCaughrean 
and Stauffer 1994) with typical mass M<j ~ 10~ 1 M Q , temperature 
~ 10 2±1 K and size ~ 40 ± 20 AU (Beckwith and Sargent 1996; see 
chapters by McCaughrean et al, Calvet et al. and Wilner et al.). The 
common existence of protostellar disks around YSOs, with properties 
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similar to those expected for the solar nebula, suggests that the condi- 
tions for planetary formation may be generally satisfied. 

Recent observational breakthroughs have led to the discovery of 
Jupiter mass (Mj ) planets around at least a few percent of nearby 
solar-type stars (see chapter by Cochran et al.). With the present 
data, we can assert that planetary formation is robust. 

In conventional planetary formation models (Pollack et al. 1996), 
the first stage of protoplanetary formation is the rapid build up of 
solid cores through the coagulation of planetesimals (Safronov 1969; 
Wctherill and Stewart 1989; Lissauer and Stewart 1993; Aarseth et al. 
1993; Kokubu and Ida 1996). When the core mass increases above a 
critical value (~ a few M®), quasi-static evolution is no longer possible 
and a rapid accretion phase begins (Mizuno 1980; Bodenhcimcr and 
Pollack 1986) leading to the formation of gaseous giant planets (Pollack 
et al. 1996; also see chapter by Wuchterl et al.). 

The details of this model are not yet fully worked out, but if we 
suppose the protoplanet can accrete gas as efficiently as possible, it 
will first take in gas in the neighborhood of its orbit assumed circular 
at radius r p , until it fills its Roche radius Rr — (g/3) 1//3 r p where q = 
Mp/M* is the protoplanet central star mass ratio, while orbiting in an 
empty annulus. The mass of the protoplanet is then given by M p = 
3 [4M d (r)/3M*] 3/2 M„, where M d (r) = 7rEr 2 is the characteristic disk 
mass within radius r, with S being the surface density. At 5.2 AU and 
£ = 200 g cm 2 , this gives a mass M p ~ 0.4 Mj. 

Further mass growth now depends on whether the disk has a kine- 
matic viscosity v capable of producing a mass accretion rate M p onto 
the protoplanet. If v is finite, then all of the mass flow through the outer 
disk should go to the protoplanet. To model the disk viscosity, we adopt 
the prescription of Shakura and Sunyaev (1973) in which v = aH 2 rt, 
where a is a dimensionlcss constant, H is the disk semi-thickness and 
n is the disk angular velocity. The most likely mechanism for provid- 
ing an effective viscosity in stellar accretion disks is MHD turbulence 
(Balbus and Hawley 1991), which produces a ^ 10~ 2 in a fully ionized 
disk. Note however that a may vary throughout the disk and be much 
smaller in its intermediate parts (see chapter by Stone et al.). 

Observationally inferred values of the disk accretion rate Md ~ 
2q-8±i jy[ Q y r -i (fjartmann e t a i. 1998) are model dependent and 
highly uncertain. Nonetheless, for such a fiducial Md, a protoplanet 
may attain a mass M p > 10 Mj within 10 6 yr. The mass of extra solar 
planets is ~ 1 Mj (see chapter by Cochran et al.). Unless these planets 
are preferentially formed in low-mass disks, their growth needs to be 
terminated/inhibited such that they are unable to accept all the mass 
that flows through the disk. In this paper we focus on disk/protoplanet 
tidal interactions as a mechanism for accomplishing this. 
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A protoplanet exerts tidal perturbations which, for M p ~ 1 Mj , 
may be adequate to induce the formation of a gap in the disk near 
its orbit and thereby start to limit accretion flow onto it (Lin and 
Papaloizou 1993). In general, the gravitational potential due to the 
protoplanet may be Fourier decomposed in the form: 



where <fi is the azimuthal angle, I and the azimuthal mode number m 
are integers and ui^ m = u + (I — m)n p /m is the pattern speed, with ui 
and k p being the angular and epicyclic frequency of the protoplanet, 
respectively. If the orbit of the planet has a small eccentricity e, then 
ipl,m tx el^ m l (Goldreich and Tremaine 1980; Shu 1984). For circular 
orbits, only I = to need be considered. 

Both outgoing and ingoing density waves are excited in the disk 
at the Lindblad resonances located at r — r L and where, for I = to, 
il = uj ± n/m (Goldreich and Tremaine 1978). Here k is the epicyclic 
frequency of the gas and is the disk angular velocity. In a Keplerian 
disk, k = ft. The ingoing/outgoing waves carry a negative/positive 
angular momentum flux measured in their direction of propagation as 
they move away from the protoplanet into the disk interior/exterior. 
The waves thus carry a positive, outward propagating, conserved an- 
gular momentum flux or wave action Fh. In most cases it is reasonable 
to assume the waves are dissipated at some location in the disk where 
their angular momentum density is deposited (see Lin and Papaloizou 
1993 and references therein). In the limit of a cold two dimensional 
disk, Goldreich and Tremaine (1978) found for a particular to that: 



The back reaction torque exerted on the disk interior (exterior) to the 
planet is — Fh (Fh)- Thus the inner disk loses angular momentum 
while the outer disk gains it. Hence the tendency to form a gap. Eval- 
uation of the total torque acting on each side of the disk requires the 
summation of contributions from all the resonances which for circular 
orbits amounts to summing over to. As to — > oo, the location of the 
resonance approaches the orbit. 

However, for a non self-gravitating disk and large to, the waves 
are sonic in character and thus can only exist further than a distance 
A ~ 2H/3 from the orbit beyond which the relative disk flow is su- 
personic. This results in a torque cut off (Goldreich and Tremaine 
1980; Artymowicz 1993a; Korycansky and Pollack 1993; Ward 1997) 
for to > r p /H. The dynamics in the coorbital region A Si 2H/3 is not 
wave-like and is considered below. Summing over over all resonances, 



= £m=0 -0i,mCOs[TO(0 - UJ^ m t)}, 
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the total angular momentum flux carried by the waves is essentially the 
same as that given by Papaloizou and Lin (1984), namely: 

H T = 0.23q 2 J:ry(r p /Hf, (3) 

which was obtained by direct calculation of the torques for a model 
disk in which the protoplanet orbited in a gap of width ~ H. 

Using equation (3) to estimate tidal torques, and then consid- 
ering the competition between viscous torques, which tend to fill a 
gap, and tidal torques, which tend to empty it, Lin and Papaloizou 
(1979a, 1980, 1986a, 1993) found a viscous condition for gap open- 
ing: q > 40^/ (flTp) . They proposed that gap formation would lead to 
the limitation of accretion but were only able to consider empty gaps. 
With the introduction of powerful numerical finite difference techniques 
and computers, recently it has become possible to study gap formation 
more fully numerically than previously (Artymowicz et al. 1998; Bry- 
den et al. 1998; Kley 1998; and see below). The first results of this 
work indicate that there is a transitional regime in viscosity for which 
a gap exists, but some accretion still occurs, but then being essentially 
switched off for small enough viscosity. 

In addition, Lin and Papaloizou (1993) pointed out that in a disk 
with very small viscosity, in order to form a gap it is necessary that 
H < Rr = (g/3) 1 / 3 r p . This thermal condition can be viewed in sev- 
eral ways. It means that the protoplanet's gravity is more important 
than pressure at a distance Rr from the protoplanet (generating large 
enough hydrostatic pressure forces would require gradients that would 
cause a violation of Rayleigh's criterion). From Korycansky and Pa- 
paloizou (1996) it is apparent that the thermal condition is required in 
order that the flow in the neighborhood of the protoplanet be nonlin- 
ear enough that shocks be produced which can provide the dissipation 
which is associated with gap opening (see below). This condition was 
found for a simple 2D-model disks with a barotropic equation of state. 

It may be possible with more complicated physics and the introduc- 
tion of three dimensional effects (Lin et al. 1990a, 1990b) to alter wave 
dissipation patterns and thus angular momentum deposition. However, 
because of the difficulty of clearing material near the protoplanet (it is 
difficult to see how linear waves could be dissipated closer than several 
vertical scale heights), it is likely that gap formation will occur much 
less readily if the thermal condition is not satisfied. 

The flow around a partially embedded protoplanet has been simu- 
lated in the low-M p limit in a shearing-sheet approximation (Korycan- 
sky and Papaloizou 1996; Miyoshi et al. 1998). These simulations 
provide useful clues on the flow pattern near the protoplanet. A global 
illustrative model is shown in Figure 1, in which we set q = 10~ 4 , 
H/r = 0.07, and a = 10~ 3 . The flow pattern can be divided into 
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three regions. (1) A protoplanetary disk is formed within Rr of the 
protoplanet with H ~ |r — r p \. In such a geometrically thick disk, the 
tidal perturbations of the host star induce a two-arm spiral shock wave 
with an open pitch angle. Angular momentum is efficiently transfered 
from the disk to the star's orbit around the planet (equivalent to the 
planet's orbit around the star) such that gas is accreted onto the core 
of the protoplanet within few orbital periods P Q = 2-k/uo. (2) An ex- 
tended arc is formed in the coorbital region near r p with gas streaming 
in horseshoe orbits around the L4 and L5 points. Because these are 
local potential maxima, viscous dissipation leads to the depletion of 
this region (Lin et al. 1987). (3) In those regions of the disk which 
have \r — r p \ > Rr, the protoplanet 's tidal perturbation results in the 
convergence of stream lines to form pronounced trailing wakes, both 
inside and outside r p . These high density ridges have been identified by 
some (Artymowicz and Lubow 1996) as gas streams along which mate- 
rial flows from the disk to the protoplanet. However, the velocity field 
viewed in a frame rotating with the companion (Figure 1) clearly indi- 
cates that gas in the post-shock region along the ridge line is moving 
away from the protoplanet. 




Figure 1. The flow pattern around an embedded protoplanet. On the left is 
the global trailing wave (the star is located toward the lower right). On the 
right is a close up view of the circulatory region immediately surrounding 
the protoplanet. Grey scale indicates surface density and the protoplanet's 
Roche lobe is shown with a dotted line. Model parameters are q = 10 -4 , 
H/r = 0.07, and a = 10" 3 . 

Simulations of gap formation need to consider model evolution over 
many orbital periods and have to deal with a large density contrast 
between the gap region and other parts of the disk. Special care is 
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needed to minimize the tendency of numerical viscosity to produce 
spurious diffusion into the gap region and thus significantly affect the 
results. In Figure 2, we illustrate the excitation and propagation of 
waves and the existence of a clear gap for a planet (with q = 1CP 3 ) 
interacting with a disk (with H/r = 0.04 and a — 10~ 3 ). 




Figure 2. Surface density distribution for a gap opening protoplanet. The 
surface density near the planet (indicated by a white sphere) has been 
reduced by four orders of magnitude. Waves are clearly seen propagating 
both inward and outward away from the protoplanet. Model parameters 
are q = 10~ 3 , H/r = 0.04, and a = 10~ 3 . 




Figure 3. The mass doubling timescale r 9 as a function of viscosity param- 
eter a for one Jupiter mass and H/r — 0.04 (left) and 10 Jupiter masses 
and H/r = 0.07 (right). As the gas viscosity is decreased, the accretion 
rate onto the protoplanet drops off sharply ultimately being determined 
by numerical effects for a < 10~ 4 . 
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In Figure 3, we plot the growth timescale M p /M p as a function of 
a and H/r for protoplanets of mass 1 and 10 Mj . We comment that 
it is not necessary for accretion to stop entirely for gap formation to 
affect the final mass of a protoplanet. All one needs is that the mass 
doubling time should become longer than either the disk lifetime or the 
time for the interior disk to accrete. In the latter case, the planet would 
approach the central star before it could increase its mass significantly 
(Ivanov et al. 1998). From Figure 3 we see that, for a ~ 10~ 3 and 
H/r = 0.04, the mass doubling time for 1 Mj is 10 6 yr. For 10 Mj , the 
same results apply for H/r = 0.07. The results in Figure 3 thus suggest 
that the tidal truncation process operates during planetary formation 
and is important in determining the final mass of the planet M p . 

II. Protoplanetary orbital evolution 

Resonant interaction between the disk and the planet can lead to 
changes in the orbital eccentricity e (Lin and Papaloizou 1993; Arty- 
mowicz 1993b). The way in which Lindblad torques cause a small 
eccentricity to grow has been reviewed by Lin and Papaloizou (1993). 
If the Lindblad torques only are considered, the disk matter is able to 
tidally pump eccentricity in a similar way to a rapidly rotating star 
(see below). However, corotation resonances in the disk (at which the 
pattern speed of a Fourier component of the tidal perturbation coro- 
tates with the disk) also need to be considered (see above). When these 
resonances exist, a torque is applied to the disk at the corotation radius 
rc where f2 = uj^ m at a rate given by: 

(4) 

r=rc 

where <r = k 2 /(2£!£) is the vortensity. The role of corotation reso- 
nances depends strongly on the ^-distribution. When X vanishes at 
some disk edge or when there is a gap near a low mass perturber, Gol- 
drcich and Tremaine (1980) showed that the action of torques from 
corotation resonances results in damping e more effectively than effects 
from Lindblad resonances could excite it. Note that this does not ap- 
ply to the situation where a protoplanet is embedded in a disk with no 
edges or gap and where tidal interaction is linear. 

Assuming dln<;/dlnr ~ 1 in the unperturbed disk, Artymowicz 
(1993b) suggested that the eccentricity of low-mass embedded proto- 
planets would be damped by tidal interaction. However, in the mini- 
mum mass solar nebula model where £ oc r~ 3 / 2 , vortensity is indepen- 
dent of r and then all corotation torques vanish (Ward 1993). Thus, 
whether low-mass embedded protoplanets suffer eccentricity damping 
(Artymowicz 1992) is model dependent. For a planet with a sufficiently 
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large mass to open a wide gap, corotation resonances would be excluded 
and e would increase on a timescale e 2 /(de 2 /dt) <~ q^ 1 (M^/Md)P (Lin 
and Papaloizou 1993) which for massive planets can be shorter than the 
inferred disk evolution timescale. 

If the growth of protoplanets is limited by the formation of a gap 
of width A s , tides may cause e to increase until e <~ A g /r p when a 
potential expansion valid for small e assumed in the above analysis 
breaks down. In a cold disk, the protoplanet's angular velocity at 
peri/ap-helion would then be greater/less than S](r p ^ A 9 ), resulting 
in a torque reversal. Such a process could provide a limit to the growth 
of e. But a protoplanet's increased radial excursion may also cause 
the gap to widen, in which case e might increase to yet larger values. 
Analytic treatments and self consistent numerical simulations in the 
large e limit have not been carried out so far. 

The characterization of extrasolar planets' orbits from forthcom- 
ing observations will provide some useful constraints on the origin of 
eccentricity. The coexistence of several planets with similar e but very 
different M p (or similar M p but very different e) is not a natural out- 
come of the excitation of eccentricity by disk tidal interaction, as this 
would tend to produce a (e, M p ) relation (Artymowicz 1992). 

Planetary orbital migration is induced by the difference between 
the inner and outer disk torques, AT, which react back on the planet, 
assumed to remain in a circular orbit, such that: 

dr, = 2AT 

dt r p flM p 1 ' 

Based on recent work indicating that AT is almost always negative, 
Ward (1997) suggested that embedded planets migrate inwards on a 
timescale ~ 1O 5 (M P /M0) _1 yr. If this process occurs in protostellar 
disks, protoplanets would migrate towards the stellar surface at r = 
i?* once M p > 1 because their growth timescale (Lissauer and 
Stewart 1993) would then become longer than this migration timescale. 
Resolution of this rapid migration delemna may require the complete 
and nonlinear analysis of the disk response to the protoplanet in the 
corotation regions which may be quite complex (see Figures 1 and 2). 
However, supposing that migration occurs and can be stopped near the 
star (see below), a scenario leading to the formation of short-period 
planets has been proposed (see chapter by Ward and Hahn). 

After gap formation, there may still be some residual accretion 
depending on the value of a (see above) . Torques are still exerted be- 
tween disk and protoplanet and the imbalance due to differing proper- 
ties of the disk on either side of the protoplanet causes orbital migration 
(Goldreich and Tremaine 1980; Lin and Papaloizou 1986b; Takeuchi 
ct al. 1996). The estimated effects of advcction of angular momen- 
tum through the accretion flow on this process are usually found to 
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be small. For small M p , the protoplanet behaves like a disk particle 
and so migrates towards the star. For larger masses, the evolution is 
slower. Nonetheless, the protoplanet is always expected to reach the 
star before it has time to double its mass (Ivanov et al. 1998). 

The above discussion naturally leads to the suggestion that short- 
period planets were formed at several AU away from their host stars, 
and subsequently migrated to their present location (Lin et al. 1996). 
However, in situ formation cannot be completely ruled out (Boden- 
hcimer 1997). 

A protoplanet 's migration may be terminated near P* if (1) the 
disk docs not extend down to P* or (2) the host star induces angu- 
lar momentum transfer to the protoplanet 's orbit via tidal effects (see 
section III). For the first possibility, interaction between the disk and 
an intense (> 1 kG) stellar magnetic field has been suggested (Konigl 
1991) as the cause for the modest observed rotation period (P* <~ 8 
days) of classical T Tauri stars (Bouvier et al. 1993; Choi and Herbst 
1996). Such a magnetic field strength may be consistent with recent 
measurements (Guenther 1997). In this scenario, the stellar field is as- 
sumed to induce a cavity out to the magnetospheric radius (r m ~ a few 
10 11 cm) where P* is equal to the local Keplerian period of the disk. 
When a protoplanet migrates interior to the 2:1 resonance of the gas 
at r m , the protoplanet is decoupled from the disk and its migration is 
temporarily stalled. The main uncertainties here are in the observed 
distribution of P* (Mathieu et al. 1998). 

III. Planet star tidal interaction 

As the planet approaches P* , the tides raised on the star or planet 
by its companion become strong enough so that their dissipation leads 
to orbital evolution. If the distribution of mass of the perturbed object 
has the same symmetry as the perturbing potential, there is no tidal 
torque resulting from the interaction. However, in general, dissipative 
processes (e.g. radiative damping or turbulent viscosity) acting on the 
tides produce a lag between the response of the star or planet and the 
perturbing potential, enabling mechanical energy to be lost and angular 
momentum to be exchanged between the rotation of the perturbed ob- 
ject and the orbital motion. Only if the system is circular, synchronous 
and coplanar (i.e. the spin of the binary components and that of the 
orbit are parallel) does the tidal torque vanish (Hut 1980). However, 
such an equilibrium state is not necessarily stable (Counselman 1973; 
Hut 1980, 1981). 

The response of a star/planet to a tidal perturbation is the sum of 
two terms, an equilibrium tide and a dynamical tide. The equilibrium 
tide (Darwin 1879) is the shape that the perturbed object would have 
if it could adjust instantly to the tidal potential. It is obtained by 
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balancing the pressure and gravity forces. The dynamical tide contains 
the oscillatory response of the star/planet. It takes into account the fact 
that gravity or g modes can be excited in the convectively stable layers 
of the star/planet and that resonances between the tidal disturbance 
and the normal modes of the star/planet can occur (Cowling 1941). 

In massive close binaries, which have a convective core and a radia- 
tive envelope, the dynamical tide cannot be neglected because tidal fric- 
tion is caused predominantly by the radiative damping of the tidally- 
excited modes (Zahn 1975, 1977; Savonije and Papaloizou 1983, 1984, 
1997; Papaloizou and Savonije 1985, 1997; Goldreich and Nicholson 
1989; Savonije ct al. 1995; Kumar et al. 1995). In that case, g modes, 
which are excited mainly near the convective core boundary, propa- 
gate out through the envelope to the atmosphere, carrying energy and 
angular momentum. They arc damped close to the surface, where the 
radiative diffusion time becomes comparable to the forcing period, thus 
enabling a net tidal torque to be exerted on the star or planet. 

In the case of solar-type binaries, which have a radiative core and 
a convective envelope, it was thought until recently that tidal friction 
could be well described by the theory of the equilibrium tide, in which 
only turbulent dissipation of the equilibrium tide in the convective en- 
velope is taken into account (Zahn 1977). However, recent studies have 
indicated that this is not the case. Even if radiative damping can be 
ignored, the torque derived using the equilibrium tide only is 4-6 times 
larger than that taking into account the dynamical tide for binary pe- 
riods of several days (Terquem ct al. 1998a). This is due to the fact 
that the theory of the equilibrium tide can in principle be applied only 
when the characteristic timescales of the perturbed object are small 
compared to the forcing period. In a solar-type star however, the con- 
vective timescale in the interior region of the convective envelope is as 
large as a month. Furthermore, because of the uncertainty over the 
magnitude of the turbulent viscosity associated with convection, it is 
not clear that the torque due to turbulent dissipation acting on the 
full tide in the convective envelope is more important than that due to 
radiative damping acting on the g modes which propagate inwards in 
the radiative core (Terquem et al. 1998a; Goodman and Dickson 1998). 
Here, we focus on the case of a planet orbiting around a solar-type star. 

Tides raised on the star by the planet can be analyzed in the limit 
that the rotational frequency of the star is usually small compared 
to the orbital frequency. Then, as a result of tidal friction, the star 
spins-up, the orbit decays (the planet spirals in) and its eccentricity, 
if any, decreases. To calculate the timescales on which this evolution 
occurs, we need to quantify the dissipative mechanisms which act on 
the tides. Recent studies (Claret and Cunha 1997; Goodman and Oh 
1997; Terquem et al. 1998a) have found that the turbulent viscosity in 
the convective envelope that is required to provide the observed circu- 
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larization rates of main-sequence solar-type binaries (Mathicu 1994) 
is at least 50 times larger than that simply estimated from mixing- 
length theory for non-rotating stars. This indicates either (1) that 
the observations are questionable, (2) that solar-type binaries are not 
circularized through turbulent viscosity acting on tidal perturbations 
(see Tassoul 1988 and Kumar and Goodman 1996 for other suggested 
tidal mechanisms), or (3) that dissipation in the convective envelope 
of solar-like stars is significantly more efficient than is currently es- 
timated (see Terquem ct al. 1998a for a more detailed discussion). 
Here we assume that circularization of solar-type binaries does occur 
through the action of turbulent viscosity on the tides, and we then 
calibrate its magnitude so as to account for the observed timescales. 
Under these circumstances, when the response of the star is not in res- 
onance with one of its global normal modes, the tides are dissipated 
more efficiently by turbulent viscosity than by radiative damping. In a 
resonance, radiative damping dominates and limit the response of the 
star at its surface. Although the planetary companion may go through 
a succession of resonances as it spirals in under the action of the tides, 
for a fixed spectrum of stellar normal modes its migration is controlled 
essentially by the non-resonant interaction. For a non-rotating star, 
the orbital decay timescale, spin-up timescale of the star and circular- 
ization timescale, in Gyr, are (Terquem et al. 1998a): 



where P Q is the orbital period. This cirularization timescale is valid 
only if the initial eccentricity is not too large. If only the convective 
envelope of the star, where tidal dissipation occurs, is spun-up during 
tidal evolution, then the spin-up timescale has to be multiplied by I c /I, 
where I c and I are the moments of inertia of the convective envelope 
and the entire star, respectively. For the Sun, I c /I ~ 0.14. 

Tides are also raised on the planet by the star. Contrary to the 
giant planets of our Solar System, Jupiter-like planets on a close orbit 
are expected to have an isothermal, and thus radiative (convectively 
stable) envelope (Guillot et al. 1996; Saumon et al. 1996). For these 
planets, both turbulent dissipation in the convective core and radiative 
damping in the envelope act on the tides. So far it is not clear which 
mechanism is more important. 
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We express the timescales associated with turbulent dissipation of 
the tides in terms of the parameter Q, the inverse of which is the effec- 
tive tidal dissipation function (MacDonald 1964). Damping of the tides 
in the convective core of the planet leads to the synchronization of the 
planet rotation with the orbital rotation on the following characteristic 
timescale, in Gyr (Goldreich and Soter 1966): 

tPp 1 (Gyr) ~ 4.4 x 1(T 13 Q 

(9) 

where R p is the planet radius, a is the semi-major axis of the orbit 
and P p is the initial value of the planet rotational period (before it 
undergoes tidal interaction with the star) . Since the rotational angular 
momentum of the planet is in general small compared to its orbital 
angular momentum, synchronization of the planet occurs before any 
significant orbital evolution can take place. Once synchronization is 
achieved, damping of the tides raised in the planet always leads to the 
decay of the orbital eccentricity on a characteristic timescale which is, 
in Gyr (Goldreich and Soter 1966): 

e<G yr) ~ 2 .8xlO-»Q^(A) 5 . (10 ) 

The tides raised on the planet do not lead to the decay of the orbit 
once the planet is synchronized. We note that Q may depend on the 
tidal frequency as seen by the planet, and therefore on the rotational 
frequency of the planet. Since t^£ c is calculated assuming synchroniza- 
tion, the value of Q in the above equation may be different from that 
used for calculating tj^ 1 . Orbital evolution of Jupiter's satellites leads 
to an estimate of Q ~ 10 5 — 10 6 for this planet (Goldreich and Soter 
1966; Lin and Papaloizou 1979b). However, it is not yet understood 
where this value of Q comes from, since turbulent viscosity arising from 
convection would produce Q ~ 5 x 10 13 (Goldreich and Nicholson 1977; 
sec Stevenson 1983 for an alternative). There is of course no reason to 
assume that the Q-value of Jupiter is similar to that of the extra-solar 
planets. First, some of these planets orbit very close to their parent 
star, and therefore are much hotter than Jupiter. Also, Q may depend 
on the magnitude and the frequency of the tidal oscillation, in which 
case it would be different for Jupiter if this planet were synchronously 
rotating on a closer orbit. Therefore we can only speculate when ap- 
plying the above formulae to extra-solar planets. 

Radiative damping of the tides in the envelope, as described above 
for massive binaries, gives rise to the synchronization and circulariza- 
tion timescales t^ 2 and i£; rc , respectively. These timescales have been 
evaluated (in particular for 51 peg) by Lubow et al. (1997). However, 
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as they pointed out, the asymptotic analysis they use is valid only if 
the initial spin rate of the planet is smaller than half that of Jupiter. 
Besides, this analysis neglects the effect of rotation on the tides, which 
is important for near-synchronous planets. 

In the context of extrasolar planets, we first comment on the mag- 
nitude of the perturbed velocity induced by the tides at the stellar 
surface. Terquem et al. (1998a, 1998b) have found that in the case 
of 51 Peg, this velocity is too small to be observed. This result is 
insensitive to the magnitude of the stellar turbulent viscosity and is 
not affected by the possibility of resonance. It also holds for the other 
extra-solar planets which have been detected so far. 

The discussion of the previous subsection indicates that the rota- 
tion of planets on a close orbit is almost certainly synchronous with the 
orbital rotation. For the star to be synchronized in less than 5 Gyr, i.e. 
t* p < 5 Gyr, a 1 (3) Jupiter mass planet would have to be on an orbit 
with a period less than 1.4 (3) day(s). We note that in the case of r Boo, 
the rotation of the star may then have been synchronized as a result 
of tidal effects. If only the convective envelope of the star is spun-up, 
these periods have to be multiplied by about a factor 2. Therefore, the 
observation of anomalous rapid rotators could give further evidence of 
the presence of close planets and provide some indication on how the 
dissipation of tides affect the rotation of the star (Marcy et al. 1997). 

The circularization timescale for the orbit is £ c i rc such that 1 / i C i rc = 
lAcirc + V*circ + V*circ- We note that, according to the expression of 
the timescales we have given above, i* rb ~ 6i* irc and t^ h > 10t£; rc , so 
that circularization can be achieved without a significant orbital decay 
taking place (in contrast to the statement by Rasio et al. 1996). So far 
it is not clear which of the timescales i* irc , t^ c and t^ c is the shortest. 
For 51 Peg, t^ c < t* irc requires Q < 2 x 10 7 (see also Rasio et al. 1996 
and Marcy et al. 1997 for estimates of t^l ). We point out that for 
a Jupiter mass planet, f* irc < 5 Gyr for P a < 3 days. Therefore, all 
the Jupiter mass planets orbiting with a period less than about 3 days 
should be on a circular orbit. This cutoff period is a lower limit, since 
i* irc is an upper limit of the circularization timescale. 

Planets found on eccentric orbits at smaller periods are indicative 
that the stellar rotation frequency (assumed zero up to now) is large 
compared to the orbital frequency Indeed, in that case tidal friction in 
the star pumps up the orbital eccentricity, opposing the effect of tidal 
friction in the planet. However, the above calculations indicate that 
the planet may have to be on a very close orbit for the tidal friction 
in the star to be able to increase significantly the orbital eccentricity. 
Whether the planet can get to such a close orbit is questionable if the 
disk in which the planet formed did not extend down to R*. 

The timescale on which the orbit decays is i* rb . The condition 
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i* rb < 5 Gyr requires P < 2 days for a Jupiter mass planet. If a 
planet has not been able to get to such a close orbit because, e.g., the 
inner parts of the disk were truncated, orbital decay may not occur as 
a result of tidal friction, and it may not plunge onto their parent star. 

However, several planets may attain short periods. After the first- 
born protoplanet emerges and migrates to a few because the for- 
mation timescale increases with r, additional, subsequently formed pro- 
toplanets may migrate inwards pushing protoplanetary cores ahead of 
their inward path until they become trapped at Lindblad resonances. 
Such resonant trapping is found for the Galilean satellites (e.g. Gol- 
dreich and Peale 1966; Lin and Papaloizou 1979b) and it occurs if the 
migration rate is sufficiently slow. It raises the intriguing possibility 
that planets temporarily parked close to the star because of absence of 
disk material could be forced to plunge into the host star. 

On the other hand, if the star is a rapid rotator (with period shorter 
than that of the orbit) and has a weak magnetic field, tidal effects trans- 
fer angular momentum from the star to the orbit with an increasing 
efficiency as the orbit is pushed out until a tidal barrier is produced 
such that the migration time and inverse orbital decay time balance. 
Tentatively using estimates derived from equation (7), an inward mi- 
gration rate of 10 6 yr could only be balanced by tidal effects acting on a 
1 Mj planet if the period was ~ 0.3 days, which increases to ~ 0.7 days 
for a 3 Mj planet and a migration time of 10 7 yr. If planets can survive 
at such short periods (more likely for massive planets and slow migra- 
tion rates) and if eccentricity growth is suppressed by tides within the 
planet, tidal effects might temporarily cause the halting of migration 
by transferring angular momentum from the star to the inner orbit and 
from there to any residual disk material via Lindblad torques and so 
on to outer planets (cf. the situation for Saturn's rings). Such pro- 
cesses have not yet been worked out in detail. However, a solar-type 
star is likely to have a small amount of angular momentum in compar- 
ison to the disk/planet system in total, and once the star slows down 
sufficiently, the inward migration must continue. 

As a planet eventually becomes engulfed by the host star, it is 
disrupted by tidal breakup, heating, and ram pressure stripping (Sand- 
quist et al. 1998). Supporting evidence for such phenomena may be 
found in the super-solar metal abundance of 51 Peg (G2V), 55 p Cnc 
(G8V), and r Boo (F7V) (Butler et al. 1997; Gonzalez 1997, 1998). 
The convective envelope for each of these stars contains a few 10~ 2 M . 
The mixing of 10-40 M of planetary heavy element material within 
the stellar envelope (Zharkov and Gudkova 1991) would lead to a signif- 
icant metal enrichment there. Since the depth of the convection zone 
of main sequence stars decreases significantly with increasing stellar 
mass, the routine engulfing of planets by their hosts might lead to a 
tendency for hotter planetary hosts to show a general overall metallicity 
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enhancement with respect to cooler ones. 

IV. Long term stability of planetary systems 

In a relatively massive disk, several giant planets may be formed 
with M p ~ 1-3 Mj and a > 1 AU. Their long term orbital stability 
determines the dynamical evolution of the system. After the depletion 
of the disk gas, mutual gravitational perturbation between the planets 
may gradually increase their eccentricities until their orbits cross each 
other on a timescale t x . Extrapolation of existing numerical results 
(e.g. Franklin ct al. 1990; Chambers et al. 1996) give t x ~ 10 12 ~ 18 yr 
for the solar system. But t x would reduce to ~ 10 8 ~ 12 or 10 5 ~ 8 yr if 
all the gas giants had M p =1 or 2 Mj , respectively. 

Since systems of massive planets formed with similar values of a 
and small e eventually suffer orbit crossing (Lin and Ida 1997), mas- 
sive eccentric planets may have acquired their orbital properties as a 
consequence of orbit crossing (Rasio and Ford 1996; Weidenschilling 
and Marzari 1996; Lin and Ida 1997). Once orbit crossing occurs, close 
encounters eventually take place. When these involve equal mass plan- 
ets, they produce a velocity perturbation V with magnitude limited 



km s 1 , where p p is the planet's internal density (Safronov 1969). This 
results in e ^ V esc /{aQ) ~ 2 (M p /M,,) 1/3 (a/1 AU) 1/2 . Thus observed 
eccentricities can be accounted for. 

Supposing that giant planets form at <~ 1 AU, massive planets 
(M p ^ 5 Mj ) with moderately high e <~ 0.3 and moderately small 
a ~ 0.5 AU can be accounted for by orbit crossing and merging. Ec- 
centric planets with a > 1 AU and M p ~ Mj can be produced by orbit 
crossing followed by ejection. Merger, resulting in e ^ 0.5, is favored 
at small r whereas ejection occurs preferentially at large r because the 
cross-section for direct collisions is independent of r while that for close 
encounters is cx V~ 4 oc (r£!)~ 4 oc r 2 (Lin and Ida 1997). Ida and Lin 
(1998) find that the orbital properties of a merged body are consis- 
tent with those of the planets in eccentric orbits around 70 Vir and 
HD 114762. 

Numerical simulations of many planet systems indicate that al- 
though some may be ejected, a residual population of eccentric planets 
(e 0.5) may remain bound to the central star at large distances (a > 
100 AU). Close encounters also excite the relative inclination up to 
~ e/2 radians. Detection of additional companions around stars with 
an eccentric planet would provide tests of these scenarios. 

Rasio and Ford (1996) suggested that the short period planets were 
dynamically scattered into a region close to the central star. Although 
the short-period planets may be circularized (see above), numerical sim- 




16 



Lin et al. 



ulations indicate that the scattering scenario would lead to a population 
of planets with high eccentricity at large distances (10-100 AU) from 
their central stars. A comprehensive search for these planets would 
provide a useful test of this scenario as against that of disk-planet 
interaction which produces orbits with smaller e.. 

V. Effects of secular perturbations due to a distant companion 

Some extrasolar planets are found in binary systems. For example, 
a planet is found around 16 Cyg B which orbits at a distance of 1.7 AU 
from the central star which is is known to have a binary companion 
16 Cyg A. It has been suggested that the high e (=0.67) of 16 Cyg B 
is excited by the gyroscopic perturbation due to 16 Cyg A (Holman 
et al. 1997). A similar effect may be caused by the gravitational per- 
turbations from other planetary bodies. Here we discuss the effect of 
secular perturbations on a long timescale such that a time average may 
be performed. The orbiting bodies may then be considered as having 
their mass distributed continuously around their orbits as in the classi- 
cal theories of Laplace and Lagrange (see Hagihara 1972 and references 
therein). Semi-major axes do not change under secular perturbation. 
However, changes to the eccentricity and inclination may be produced. 

If orbits arc coplanar, in general only modest changes to eccen- 
tricities are produced if the perturbing bodies are widely separated. 
An exception to this might occur if secular resonances sweep through 
the system because of a changing gravitational potential due to disk 
dispersal or changing stellar oblateness (eg. Ward 1981). However, if 
orbits are allowed to have high inclination, large eccentricity changes 
can occur to the orbit of an inner planet perturbed by a distant body 
or binary companion. 

To consider this gyroscopic perturbation effect, let us consider the 
simplest case of the motion of an inner planet with mass M p perturbed 
by an outer companion with mass M' assumed to be in a circular orbit 
with semi- major axis a'. The outer companion is assumed to contain 
sufficient angular momentum that its orbit remains fixed defining a 
reference plane to which the inner planet orbit has inclination i. We 
note uj a the longitude of the apsidal line measured in the orbital plane 
from the line of nodes. 

If the distance of M p from M* is r, its speed is v, the angle be- 
tween its position vector and that of M' is O, and only the dominant 
quadrupole term in the interaction potential multipole expansion is 
retained, its motion is governed by the Hamiltonian (Kozai 1962; Hag- 
ihara 1972; Innanen et. al. 1997): 




GM*M, 



r 



p 



GM'M p r 2 
2^3 



(3cos 2 9-l). (11) 
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After performing a time average, appropriate for secular perturbations: 

GM*M P GM'M„a 2 r - 2 „ 2x 

= 2a ^[(Scos^-l^ + Se 2 ) 

+ 15e 2 sin 2 icos2w a ] . (12) 

The fact that both H and the component of angular momentum parallel 
to the outer orbital axis, Lcosi = J GM,o(l — e 2 ) cos i, are constant 
enables elimination of i and u> a , so that a complete solution can be 
found from the equation for the rate of change of e, which can be 
derived from the canonical equation dL/dt = —dH/duJa, in the form: 

de 15nM'a 3 esin 2 i sin2u; a 

~dt = 8M^3 ' { ' 

with 27t/n being the orbital period of the planet M p . 

The above equation indicates that e can oscillate between extremes 
occurring when sin2w a = 0. There has been particular interest in find- 
ing conditions under which e can start from very small values and then 
increase unstably to values & 1. To examine conditions for this to 
occur, we use H=const and Lcosz=const to find u a and i in terms of 
e. Assuming initial values i — io and e = eo, when uj a = u>o, for in- 
finitesimally small but non zero values of eo that can be neglected, we 
find: 

1 — 5 cos 2 i — e 2 

cos2w a = — 2~ T- 14 

5(1 - e z — cos 2 1 ) 

From the above, we see that when e — > 0: 

(l-5cos 2 i ) 
COS2 " a= 5(l-cos^ o) - (15) 

Clearly, for solutions of the type we seek, we require cos2uu a > —1, 
which requires that the initial inclination exceed a critical value (see 
Innanen et. al. 1997) such that: 

cos 2 i <3/5. (16) 

If this inequality is satisfied (if it is not then eo cannot be neglected), 
equation (13) indicates that e grows from a very small value up until 
the value obtained by setting ui a = ir/2, namely: 

e 2 = 1 - 5cos 2 i /3, (17) 



attained when cos 2 i = 3/5. A range of eccentricities may be generated 
in this way. However, values close to unity may be attained for initial 
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inclinations close to 90°. There are then solutions which have large 
amplitude oscillations in eccentricity (Kozai effect) and adding back 
in the non secular terms can lead to chaotic behavior. The large ec- 
centricity changes occur independently of the size of the perturbation. 
However, the characteristic timescale for the eccentricity changes to 
occur is given from equation (13) as ~ (M*a l3 )/(M'a 3 ) orbital periods 
of the inner planet, which is long for small perturbations. 

But note that the effect requires the motion of the apsidal angle 
uj a be governed only by the perturbation considered. Other effects may 
disrupt this, such as general relativistic corrections (e.g. Holman et al. 
1997) and the oblateness of the central star. Both could be considered 
in the discussion below, but here we limit ourselves to considering the 
effects of relativistic apsidal precession which may lead to constraints 
on the orbital elements of planets with short periods. 

To incorporate relativistic apsidal precession to lowest order, we 
modify the Newtonian potential such that keplerian elliptic orbits are 
induced to precess at the correct rate. We modify the potential due to 
the central mass, such that: 

-GM* -GM„ f i + 3GMA ^ 



with c being the speed of light. The additional term can be added into 
the Hamiltonian H and averaged so that equation (12) becomes: 

H^H- i{GM ^ . (19) 
2a 2 c 2 v /(l -e 2 ) 

The same procedure as described above can be used to find co a and i 
in terms of e and hence conditions that must be satisfied in order that 
large values of e might be attained. The condition analogous to (16) 
for the Kozai effect to work is: 

cos 2 i < (3-2/)/5. (20) 

Here the parameter: 

'-mm 

measures the importance of relativistic precession relative to the per- 
turbation due to the outer companion. From this, if / > 3/2 the effect 
is completely suppressed. When the Kozai effect occurs, the maxi- 
mum possible eccentricity that can be generated is obtained by setting 
uj a = tt/2 and io — 90°. This gives: 



e 2 = 1/2 + V(4//3 + l/4) - 4//3, (22) 
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which indicates that, for / = 0.1, e as high as 0.99 can be generated 
for i = 90°. The parameter / may be conveniently expressed as: 

Here, P Q and P^ denote the orbital periods of M p and M' respectively 
and q' = M'/M*. 

Planets with periods less than around 3-4 days are likely to have 
their orbits circularized as a result of tidal effects (see above). From 
the requirement / > 3/2, we see that for inner planets with somewhat 
larger periods, but still on the order of days, binary companions with 
mass ratio of order unity need orbital periods shorter than about 40 yr 
in order to pump significant eccentricity. If a massive planet with 
q' ~ 10~ 2 is considered, it must orbit with a period ten times shorter 
<~ 4 yr. Such an object should be readily detectable. 

We comment that for a given distant binary companion with incli- 
nation close to 90° to the planet orbit, and no other perturbers, there 
is a planet with a period long enough such that the Kozai effect works 
such that it could produce eccentricities close enough to unity that the 
planet has arbitrarily close approaches to the central star. Then tidal 
effects may act towards circularization so decreasing the distance to 
apoapse until a circular orbit is produced. The final period would have 
to be a few days. Alternatively, should the closest approach to the star 
be beyond the largest radius at which orbits can be circularized, the 
orbit would retain a Kozai cycle with high maximum eccentricity. 



VI. Summary 

Here we summarize the various processes discussed in the paper, 
their influence on the evolution of planetary systems, and their obser- 
vational implications. 

a) Planet-disk tidal interaction excites potentially observable spiral 
density waves and can create clean gaps in the disk, limits mass growth 
and drives orbital migration. Growth limitation by gap formation to- 
gether with inward migration may naturally account for the upper limit 
in the observed mass ratio distribution of extra solar planets provided 
they were formed in disks with H ^ 0.1r p and a & 10~ 3 . Since the 
critical mass for gap formation is an order of magnitude larger than for 
dynamical gas accretion, a bimodal mass distribution with a depres- 
sion between these masses may be anticipated. After gap formation, 
massive protoplanets are expected to migrate with the viscous evolu- 
tion of the disk until they encounter a stellar magnetospheric cavity 
or tidal barrier. The existence of short period planets in CTTS would 
imply that the local disk viscous evolution timescale is shorter than the 
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typical age (~ 10 6 yr.) Orbital migration may lead to stellar consump- 
tion of protoplanets with surviving planets locked into commensurable 
orbits. Contamination of their host stars (during the latter's main se- 
quence evolution) by plunging planets may have led to a relatively high 
metallicity. Planet-disk interaction is likely to excite only small e ^ 0.3 
and lead to a relation between M p and e which could be observed. 

Stellar tides raised on short-period planets should rapidly synchro- 
nizes their spin. Short period planets with M p ~ Mj and P & 3 — 4 days 
(e.g. r Boo) should have their orbits circularized by tides raised on a 
slowly rotating star within a few Gyr. In the case of massive plan- 
ets, the stellar rotation may also be synchronized with the orbit. But, 
only very short period planets (P & 2 days) are expected to undergo 
significant orbital decay. 

c) Long term gravitational interaction between planets with initially 
small e leads to orbit crossing on a time scale which is sensitively de- 
termined by their masses and initial separations. Subsequent close 
encounters between comparable masses can lead to high orbital eccen- 
tricity and planets scattered into extended orbits. This mechanism 
may provide a supply of planets with a range of short periods, some of 
which may undergo tidal circularization and/or plunge into the star, 
together with outlying accomplices. The latter could be imaged by the 
next generation of IR interferometers. 

d) Secular perturbations by distant binary, or massive outlying 
planetary, companions in relatively highly inclined orbits can also excite 
high e through the Kozai mechanism. Potentially this process might 
produce a supply of short period planets (P ^ 2 — 3 days) that have 
undergone tidal circularization. But because of the very high relative 
inclinations required, it is likely to work only in a small number of 
cases. 

This work is supported by NSF and NASA through grants AST- 
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